Potential Errors in a Scheme of Universal Quantum Gates in Kane's Model 
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We re-investigate a plausible proposal for universal quantum gates in Kane's model, in which the 
authors assumed that electron spin is always downward under a background magnetic field and the 
value of controlling parameters is varied instantaneously. We demonstrate that a considerable error 
appears, for example, in the X rotation. As result, the controlled operations don't work. Such a 
failure is caused by improper choice of the computational bases; actually, the electron spin is not 
always downward over time during quantum operations. 



PACS numbers: 03.67.Lx 



INTRODUCTION 



Various physical models for quantum computation 
have been proposed. Among them, a nuclear magnetic 
resonance(NMR) quantum computation by liquid state 
NMR 0, |2 is the most successful, because several 
quantum algorithms have experimentally worked, using 
five or seven qubits (see Ref. jij for additional references) . 
On the other hand, Gershenfeld and Chuang Q pointed 
out that it was not possible to achieve beyond ten qubits 
using it. 

Among many proposals for realizing quantum com- 
putation with many qubits El, j Kane [j| proposed 
quantum computing by nuclear spins in a semiconduc- 
tor. In this proposal, the nuclear spin of a dopant atom, 
31 P, implanted into a silicon substrate is a single qubit, 
and quantum states are measured by detecting currents 
of spin-polarized electrons around 31 P. Although the 
implantation of the atoms with requisite accuracy and 
the detection of the single electronic charge motion are 
very difficult even in the present state of the technol- 
ogy, Kane's proposal could achieve quantum computa- 
tion by the use of many qubits, applying the several 
existing microfabrication techniques of semiconductors. 
Moreover, several important experimental techniques, for 
example, the single-ion implantation method Q and the 
single-electron transistor 0, 0, El , have been devel- 
oped steadily. Hence, it is necessary to re-investigate the 
feasible conditions required by Kane's original proposal, 
based on detailed theoretical analyses and the present 
state of experimental techniques. 

Several investigations of the schemes of implement- 
ing quantum gates in Kane's model have been reported. 
Wellard et al. [12] numerically derived an effective Hamil- 
tonian and proposed a nonadiabatic controlling scheme 
for a controlled-NOT (cnot) gate. It should be noted that 
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the quantum computational bases chosen are not eigen- 
states of their Hamiltonian. Fowler et al. 13] showed an 
adiabatic controlling scheme for a cnot gate and stud- 
ied the errors under the presence of dephasing. Their 
proposal is based on Goan and Milburn's study, which 
is explained briefly in Ref. 0] • Hill and Goan |15| pro- 
posed a nonadiabatic scheme for arbitrary single qubit 
operations and controlled operations (e.g., a cnot gate, 
a swap gate, and a controllcd-Z gate), and studied the 
effect of dephasing. 

In particular, Hill and Goan's results are very impor- 
tant, as they make it possible to construct a set of uni- 
versal gates for quantum computation |16| . They use the 
assumption that spin of the donor electrons is always 
downward, even if quantum gates are carried out. Fur- 
thermore, a nonadiabatic, or instantaneous, controlling 
process of parameters is an ideal one from both theo- 
retical and experimental points of view. Therefore, it is 
necessary to discuss the validity of such assumption and 
idealization. 

In this paper, we show that potential errors should ex- 
ist in Hill and Goan's schemes for the gates including 
spin-flip operations (e.g., X rotations and cnot gates); 
there is a large difference between their aimed state 
\4>target) obtained by quantum operations and the phys- 
ical state \<j)phys) obtained by a more realistic time evo- 
lution. We think that they didn't properly take account 
of the hyperfine interaction (HY) between electron spins 
and nuclear spins in their composition of the qubits sys- 
tem. 

This paper is organized as follows. First, we review 
Kane's model and, taking the strengths of all interac- 
tions in the model into account, explain how to choose 
our computational bases in Sec.[H] After the completion 
of choosing them, we discuss the schemes of several im- 
portant quantum gates, phase-shift operations, spin-flip 
operations, and controlled-Z operations in Sec. lIIII The 
main results are shown in Scc. lIIIBl Finally, we summa- 
rize our results in Sec.lIVI 
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II. KANE'S MODEL 

Kane's model for two qubits is described by the Hamil- 
tonian 

2 2 
H(t) = J2 H\t) + J(t)a le • <x 2e + ]T H Ut) > (1) 

i=l i=l 

where H l (t) = ~g n p n Ba l z n + n B Ba\ e + A t {t)<r te ■ cr ln , 
H l ac (t) = B ac m(t) ■ (-g n fi n cr in + u B cr le ) and m(t) = 
(cos(<J oc t), — sin(uj ac t), 0). The Pauli matrix is written 
by Cfc (k — x, y, z), and the superscript ie (in) in it in- 
dicates the electronic (nuclear) spin of the dopant atom 
implanted in the i th site. According to Ref. 5], the value 
of the uniform static magnetic field in the direction of the 
z axis is 2.0 T (i.e., B = 2.0 T). Then, the value of the 
Zeeman splitting energies for the electron, p B B, and the 
nucleus, g n /j, n B, are 0.116 meV and 0.071 x 10 _3 meV, 
respectively. The HY on the i th site between the nu- 
clear spin and the electronic spin, Ai(t)cr le ■ <j %n is locally 
changed by controlling the voltage in the i th A-gate lo- 
cated right above the i th dopant atom p| . The typi- 
cal value of the HY is 0.121 x 10~ 3 meV (= A ) when 
the voltage in the A-gate is vanishing, and the mag- 
nitude of Ai decreases as it increases p . Notice that 
the nuclear Zeemann splitting energy and the HY are 
almost the same order of magnitude. This causes a non- 
negligible mixing between the electron spin down and up 
in the eigenstates of the Hamiltonian (see below). The 
electronic exchange interaction (EE) between the 1st and 
the 2nd electronic spins, J(t)cr 1 e • er 2e , is locally changed 
by controlling the voltage in the J-gate located between 
the neighboring dopant atoms Q . We assume that it 
will be possible to change the magnitude of the EE from 
J = to J ~ hbB/2 The third term in Eq. Q is 
a time-dependent interaction between the magnetic field 
perpendicular to the z axis and spins. The typical value 
of B ac is assumed to be B ac ~ 2.5 mT [5|,[l5|. We write 
the eigenstates for a\ n as |0), and |1) 4 (a™\Q)i = |0), 
and a" l \l)i = — |l)i), and those for <j\ e as and \{)i 
(<|T)i = |T)iand <%\l)i = 

To investigate the time evolution of the system in de- 
tail, we first choose the states to represent a qubit from 
the eigenstates of H 1 (t ) with a fixed value of Ai (t ) . These 
eigenstates compose the temporal bases in the adiabatic 
calculation. We abbreviate the argument t in the equa- 
tions of its eigenvalues and eigenstates (we will recover it 
if necessary). Introducing e = fi B B + g n fi n B, the eigen- 
values El(Ai) (k — 0, 1, 2, 3 from the below) are given 
as follows: 

E'o(Ai) = -A t - yJe* + 4Al 
E\(Ai) = -u B B + g n u n B + A u 
Ei(AA = -A i + y /e* + 4Al 
Ei(Ai) = fi B B - g n fi n B + Ai. 
Defining cos(9 4 = (e + ^ e 2 + 4A 2 ) / N u sin (9, = 2A i jN i , 



and Nf = 2y/e 2 + AA 2 (e+^e 2 + 4A 2 ), their correspond- 
ing eigenstates are also given by 

K(^)>* = -sintfi |T><|1>< + oob«* U>»|0> 4 , 
\u 2 (Ai))i = cos ^|T)i|l)i + Biiifl i U) i |0) i> 

l«s(^))i = |T)i|0)i. 

The operator S l z = [a\ e 4- a 1 ™)/ 2 commutes with 
H\ Si\uo(Ai))i = Si|u2(^))< = 0, S*MA)>i = 
-\ui{Ai))i, and Sl\u 3 (Ai))i = |u 3 (^;))i- The ground 
state is always \uo(Ai))i. The energy differences A l k0 = 
E l k (Ai) — Eq(Ai) between the k th excited states and 
the ground state are estimated as follows: A^ ~ 2Ai + 
2g n n n B, A| ~ 2e, and A\ ~ 2e + 2A t - 2g n p n B. The 
values of two terms in A^ are of the same order of mag- 
nitude as HY; it is possible to control the phase differ- 
ence between |iio(j4j)), and \u\(Ai))i, varying the mag- 
nitude of the HY, Ai(t). Note that these two states are 
robust against this controlling process because they are 
the eigenstates of temporal Hamiltonian H l (t). Conse- 
quently, we choose the quantum computational bases for 
a single qubit as |uo(Aj))$ and \ux(Ai))i. In particular, 
we write jO)^ = \u (A )) l and = \ui(A ))i when 

Ai(t) = Aq. 

Next, we discuss polarization of the electron spins. 
Taking account of smOi ~ 10 -3 for Ai = Aq, we may 
regard \uo(Ai))i ~ |J.)i|0)j, and the electrons are almost 
polarized in the downward direction of the z axis. If the 
system is a thermal equilibrium state with the tempera- 
ture T = 100 mK (/T 1 = k B T ~ 8.62 x 10" 6 eV), which 
is a typical value in Kane's model [5j , the ratio of the 
number of upper polarized electrons in the i th site, n* e , 
to the number of downward polarized ones in it, n 1 ^ , is 
estimated as follows: 

nf _ Tr(|t). t (tM g ') 
nf ~ TSr(U)iU|e-^) 

e -^A*3o + cos 2 6i e-P^o + sin 2 ^ 
sin 2 9i e~P A io + e ~' 3Al io + cos 2 6i 
~ 10~ 5 . 

Here, we have used the following estimation: sin0j ~ 
10" 3 , 2(3t ~ 2fi B B/k B T ~ 10, and 2/L4 ~ 10" 2 . The 
result, nl e /n l f ~ 10~ 5 , also suggests that the electron 
spins should be |J,)j. However, the above consideration 
is based on the static physical property; it is nontrivial 
that it could sustain the downward polarization through 
the temporal controlling process. 

Now, we address the preparation of an initial state. 
The thermal equilibrium state of the system p eq = 

®i(e~ /3HZ /Tr e~^ HZ ) is not a pure state, even if the 
temperature T is 100 mK. In fact, we find p l eq ~ 
0.508 \uo(A i )) i {u {A i )\+ 0.492 for T = 

100 mK, A t = A Q , J = 0, and B = 2.0 T. The initial- 
ization problem might be solved by the standard scheme 
in liquid NMR (e.g., by the use of effective pure states). 
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Throughout this paper, we investigate only whether a 
given pure state evolves to a desired state by controlling- 
processes. Such consideration is meaningful, because any 
arbitrary mixed state is represented by suitable linear 
combinations of pure states, and the control of a quantum 
system according to our expectations will be necessary to 
overcome initialization problems. 
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FIG. 1: Configurations of the He energy level around the 
level crossing point (J ~ (j,bB/2). Four solid lines repre- 



sent the lowest eigenvalues of Hq' x \ H, 



r(- 



H, 



and 



H, 



(-2,1) 



respectively, which are associated with the compu- 



tational bases. Two dashed lines represent the second lowest 
eigenvalues for each Hg' 1 ' and Hq 1 ' 1 '. 



Hence, we calculate analytically the diagonal form of 
Eq. J2J. We show the explicit forms of its eigenvalues 
and eigenstates in AppendixIXI We explain the rela- 
tionship between the eigenstates of He if) and the com- 
putational bases. If the value of Jit) is vanishing, we 
hnd that = |0) L , 1 |0) L , 2 , |«+) = (|0) L , 1 |1) L , 2 + 

|1>l,i|0)l, 2 )/\/2, = (|0) L , 1 |1) L , 2 - |1)l,i|0)l, 2 )/V2, 
and \vf) = |l)i,i|l)i,2 are the eigenstates of fQ^ (t), 



H, 



and H c 2 ' +1 \t), respectively. 



Therefore, these four states play an essential role in the 
scheme of controlled operations. In Fig.Q] we show that 
the eigenvalues of submatrices related to the computa- 
tional bases around the level crossing (J ~ hbB/2), 
varying the magnitude of the EE, J. The eigenvalues 
of H^' p ' are specified by {k = 1, 2 for (s, p) — 

(0, 1), (-1, 1), and k = 1 for (s, p) = (-1, 1), (-2, 1)). 



The energy difference E\ 



(-1,1) 



(-i,-i) 



tive for J < /ibB/2. As J increases, the levels E\ ' and 



is tiny and posi- 

(0,1) 



E 



(-1,-1) 



approach E^' 1 and E$ , respectively. 

III. THE QUANTUM GATES 
A. The phase— shift operations 



Finally, let us consider the eigenvalue problem of the 
two qubit Hamiltonian 



Hdt) = J2 Hi (t) + J®* 1 



r 2 e 



i=l 



fixing the values of A\{t), A 2 if), and Jit). The Hamil- 
tonian He it) is composed of the first and the second 
terms in Eq. and plays a central role in the controlled 
operations. Hereafter, we write \iji)i\<j>)i as \ip4>)i (e.g., 
= || l)i). We introduce the total spin opera- 



tor in the direction of the z axis SI 



ial 



a 



le 



o\ n + <r 2e )/2 and the parity operator, P defined as fol- 
lows: P*a l a a P = <j\ a and P ] a 2 a a P = a\ a (a = e, n 
and a — x, y, z). The operator P is an idempotent op- 
erator: P 2 = 1. We find that Sl ot always commutes 
with both He it) and P, and P commutes with He it) if 
At(t) = A 2 (t). As result, in the case of Aft) = A 2 (t), 
the Hamiltonian He it) becomes a block diagonal ma- 
trix whose blocks are characterized by the eigenvalues 
of Sl ot and P, s (e <r(S* ot ) = {2, 1, 0, -1, -2}) and p 
(e c(P) = {1, —1}), respectively: 



H c it)=®H^\t) 



(2) 



s, p 



Noh> thai, Si ot H { c s,p \t) = 



H^ p \t)Sl ot = sH^ p >(t) 

and PH%' p \t) = H ( f' p \t)P = pH { f' p \t). A simi- 
lar result is given by Berman et al. |17| , but our com- 
putational bases are different from theirs. The sub- 

r(s>p) 



(«.p)/ 



We achieve the phase-shift operation for the i th qubit, 
varying Aft) adiabatically with the other parameters 
fixed; Jf) = 0, B ac = 0, and Ajf) = A (j f i). 

First, we calculate the time-evolution operator for the 
i th qubit: 



Uft z ) = T jexp 



H l f)dt 



(3) 



The symbol T means the time-ordered product. We 
define the operation time for the phase-shift operation 
as tz{> 0). Introducing a dimensionless parameter a 
(0 < a < 1) and a dimensionless variable r = t/tz, 
we assume the time dependence of Aft) as follows (see, 
Fig.©: 



Aft) 



AofAf) 

A f A (l-r) 



< r < i 



(4) 



The function /(r) in Eq. Q is defined as 



/a(t) = 



1 - 8ar 2 



-eir-- 



1-a- 



8a(r-- 



and O(r) is the Heaviside function. The initial state for 
the i th qubit will be spanned by |0) L,i and 1 1) L,i, because 
AfO) = A . We also find that \uo(Ai(tz)))i = \0)l,% and 



\ufAftz)))i 



because we assume that the time 



matrices Hq ft) in Eq. (O are 4x4 ones at most. dependence of Aft) is given by Eq. Q. Consequently, 
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FIG. 2: Profile of Ai(t) (a = 0.6). 

Eq. J2J), due to the adiabatic theorem [lil Hs| . is written 
by 

I/ i (tz)=e-« , «|0)i, i (0|+e-^|l>x,i<l|+i2 i J (5) 



where Sq = tzh 1 J Q — Ai(s) 



ds and 



The third 



g„fi n B - HbB + Jg 1 Ai(s) ds 

term in Eq. (J5J, R l includes \u2(Ai(tz)))i(u (Ai(tz))\ 
and \u (Ai(tz))}i(u2(Ai(tz))\ mediated by the level 
crossing effect in the course of operation. By the defi- 
nition of a matrix norm 20] for an n x n complex matrix 
M = (M ab ), \\M\\ ee maxi< a , h <„ |M &|, is estimated 
as the order of h/2etz, and later on we find that it's very 
small. A criterion for the validity of the adiabatic ap- 



proximation is given as follows \li 



max 

0<T<1 



(t4(r)|dH*(r)/dr|tij(r)) 



E\{t)-EI{t) 



< 1. 



(6) 



In this case, we have a condition, 



eA 



AA 2 



a < 1. 



The inequality © is always fulfilled, because 0(2A^/e) ~ 
10" 3 and O(a) ~ 1. 

Next, we calculate the time-evolution operator for the 
remaining j th qubit (j ^ i). Taking account of the 
constancy of Aj(t) = Aq and the initial states for the j 
th qubit spanned by \0}L,j and |l)z,j, we simply obtain 
the following result: 



W(t z ) = e-<\0) L , j {0\ + e-^\l) L , j {l 



(7) 



where Sq = t z h 1 (— A — yfe 
tz^ 1 



4A 2 ) and §{ = 

{gnUnB - H B B + Aq). 

Now, we construct the phase-shift operation (i.e., a Z 
rotation) with the given angle 9z for the i th qubit. In 
the first place, the phase difference between \0)l,i and 
|1) la should be equal to 9z\ Sq — 6\ — Qz + 2mr (n e Z), 
due to Eq. iJSJ. Secondly, for all qubits except for the i th 
qubit, no phase differences have to exist; Sq — S± — 2m/n 
(m 6Z, 0, and j ^ i), due to Eq. Q. Using these 
two conditions, we obtain the following equations: 



J 



e - 2g n ^ n B + J 1 \-2A t {s) - ^TlMsf ds 



2nn 



tz 



e - 2g n (j,„B - 2A - ^e 2 + AA^ 
2nmh 



2w/k 



e - 2g nf i n B - 2A - ^ e 2 + AA 2 

I 



(8) 
(9) 



We find that Eq. (jHJ includes only one unknown parame- 
ter a, except for the free parameters m and n. We should 
choose suitable values for them so that the solution of 
0(a) ~ 1 exists for any 9z; actually, we choose m = — 5 
and n = —6. We numerically solve Eq. JSJ by the bi- 
section method. Then, we get the operation time tz by 
substituting the value of m for Eq. © ; its value is inde- 
pendent of the value of Qz- 

We show the results in TableQ] The minimum value of 
Ai(t) is given by Ao(l — a). The phase-shift operations of 
the angle 9z = tt/4 (i.e., the n/8 gate) and n/2 (i.e., the 
phase gate) are performed by reducing the HY by up to 
about 66% during the operation time, which is about two 
times longer than that of Ref. 0] . This difference comes 
from our chosen profile function of Ai(t), valid for the 
adiabatic approximation. The value of h/2etz is approx- 



TABLE I: 


Calculated values of a and tz with m 


= —5 and 


n = —6. 






0z 


a 


tz (/is) 


tt/4 


0.598 


0.05 


7r/2 


0.664 


0.05 



imately estimated as Ao/4\m\ire; Ao/4\m\ire ~ 10~ 5 for 
to = —5; the correction to the adiabatic approximation is 
very small. The choice of m and n is not unique. In fact, 
we find that other solutions of Eq. © exist if m < —5 and 
n < —6. The operation time tz for (to, n) = (—5, —6) 
is the shortest of the solutions. We confirm also the re- 
alization of the desired phase-shift gates by solving nu- 
merically the time-dependent Schrodinger equation for 
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the resultant parameters. 

We summarize the results for the phase-shift opera- 
tion. We have achieved it by adiabatically controlling the 
magnitude of the HY, Ai(t), with J(t) = 0, Aj(t) = A a 
(j =/= i), and B ac — 0. Its errors are very small, as the adi- 
abatic approximation works well. The phase-shift oper- 
ations relevant to universal quantum gates, the 7r/8 gate 
(0z = tt/4) and the phase gate (9z = tt/2), are realized 
by reducing the value of Ai(t) from A to about 0AA o at 
most during the operation time tz(— 0.05 /is). Accord- 
ing to , it might be possible to decrease the magnitude 
of the HY by up to 50%. This means that further devel- 
opment of the experimental techniques might be needed 
to fully realize them by our scheme. 



B. The spin— flip operations 



where r = t/t' x , tx = tx/t' x , and r sp — 1 + Tx- The 
parameter a in /a(t) is given by a = 1 — A/Aq. In 
addition, the transverse time-dependent magnetic field 
is globally applied, while we keep Ai(t) = A (the second 
step). 

A undesired phase difference between |0)i and 
is caused by the adiabatically varying processes during 
Ai(t) (the first and third steps). This phase difference 
can be canceled by carrying out a suitable Z rotation, 
as described in Sec. lIII Al before and after the first and 
third steps respectively. 

Let us discuss the Hamiltonian in the second step. Af- 
ter instantaneously turning on the transverse magnetic 
field, the Schrodinger equation for the i th qubit is writ- 
ten from Eq. (JTJ): 



Let us consider the scheme of the spin-flip operation. 
We reproduce the spin-flip operation in Ref. [15J by ap- 
plying the transverse time-dependent magnetic field (see 
the third term in Eq. Q), using an approximation. The 
idea is based on a standard technique in NMR 0, . 
Here, we assume that the transverse magnetic field is 
turned on or off instantaneously. We show that a serious 
error exists in Ref. |Tf| , even if such an ideal control- 
ling process is realized. Hereafter, we concentrate on the 
scheme of an X rotation for the i th qubit. 

First of all, we explain the controlling processes of the 
EE, the HY, and the transverse time-dependent mag- 
netic field. We keep J(t) = and Aj(t) = A (j ^ i). 
By contrast, we adiabatically change the magnitude of 
Ai(t); in the first place, we decrease the value of Ai(t) 
up to a suitable constant value A(< Aq) (we explain how 
to choose the value of A below) during the time interval 
t'x/2 (the first step), keep Ai(t) — A during the time in- 
terval tx (the second step), and finally increase it from 
A to Aq (the third step). Thus, the time dependence of 
Ai(i) is given by 



Mt) = 



A f A (r) 
A 

AafA(T sp 



< T < 
1 
2 



I) 



T ) { T sp- 



< t < r. 
l 

2 



(10) 



(11) 



We solve Eq. in the rotating frame |2ll |23 | ; we 
remove the time dependence of the Hamiltonian in 
Eq. Ijlll) by applying the unitary operator D z {t]UJ ac ) = 
exp [~ioj ac t(al e +a l z n )/2]. Introducing \ip r ot(t))i = 
D z (t;uj ac )\ilj(t)}i, we obtain the following Schrodinger 
equation in the rotating frame from Eq. (|llfl : 



ih—\i> rot (t))i = /f; ot |Vvert(*))i, (is) 

at 



where H l rot = fiw e ^ e /2 - hw n (j z n /2 + Aa ie ■ er ln + 
\i B B ac o % x e - g n ii n B ac (j l x n , hu e = 2fi B B + hw ac , and 
huj n = 2g n /i n B — huj ac . We naturally expand \ip r ot(t)}i 
and H % rot by {\uk(A))i}. Let us introduce the matrix 
representation for an operator M as follows: M a b = 
l {u a {A)\M\u b {A)) i (a, b = 0,l, 2, 3). Thus, the Hamil- 
tonian in the rotating frame, B l rotl is written down as the 
sum of a block diagonal matrix H l d and a off-diagonal 
matrix i/L„ ■ 



< T < 



spj 



H ro t — H d + H m i x 

( Eq(A) -u e B ac 
-v Bac E\(A) - hw ac 

\ 






E%(A) -v- e B ac 
-V-gB ac El(A) + hw ac J 



[i-eBac 

fi e B ac 

n e B ac 

li- e B ac 



(13) 



Here, we define fie = ^bcos9 — g n fJ- n sin 9 and vq = 
(iB sin.0 + g n [i n cos 9, and the mixing angle 9 is defined 
for 9i when Ai it) = A. Whenever one can safely assume 
that the off-diagonal Hamiltonian, B % mix , is negligible, 



the computational bases, \ua(A))i and \ui(A))i, do not 
mix with the remaining irrelevant states, \u2(A))i and 
\u3(A))i, because of the block diagonal Hamiltonian, H l d . 
In most studies up to now, this approximation has been 
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taken for granted. 

However, it is nontrivial whether we should omit H l mix 
in solving Eq. 1)12(1 . Notice that fig is larger than vg, 
unless is close to ir/2. We show that a considerable 
error is caused by neglecting H l mix below. 

Before discussing such an error, we reproduce the re- 
sults in Ref. |Tfij by means of our computational bases. 
We solve Eq. JT^J, omitting H % mix . We choose the value 
of A so as to satisfy the Larmor resonance condition cor- 
responding to the computational bases, 



Tw ac = E{{A)-Ei{A). 



(14) 



Taking account of Eq. i|14|) and the initial state \<p)i ex- 
panded by \u (A))i and \ui(A))i (i.e., the state just be- 
fore turning on the transverse magnetic field) , we obtain 
the resultant state at tx in the rotating frame: 



\iprot(tx))i 



p-itxE^A) e i0 x a x (A) I 



(15) 



We define a x {A) = |t*o(A)>i<tti(A)| + \ui(A))i(uo(A)\ 
and Ox = veB ac tx /h. Accordingly, 

m x )) l = e-^ E o(A) D t z ^ actx) ^e x a x (A)^ )t _ (16) 

Under the Larmor resonance condition (|14|l for the i th 
qubit, the time-evolution of the remaining j th qubit 
(j i) is approximately expressed by a Z rotation. More- 
over, both the overall phase and D\(uj ac tx) in Ea. i|ltj|) 
are canceled by carrying out suitable Z rotations. Con- 
sequently, we make only the state of the target qubit flip 
due to Eq. l(TH|l . In this way, we have reproduced a result 
corresponding to that achieved in Ref. [l5j ■ 

Now, we demonstrate the existence of a considerable 
error in the above scheme. We solve Eq. 112|) without 
omitting H % mix . After this calculation, wc compare the 
exact solution with the approximate solution (|15fl by cal- 
culating the following quantity: 



HI 



(17) 



This is just the fidelity [lH] between the aimed state 
e~Ti Hdtx \4>)i obtained by an X rotation and the phys- 
ical state e~i H --t tx \(t))i obtained by Eq. (H2). If the 
value of F(\(j)}i) is adequately close to one, the dynam- 
ics generated only by H l d is a good approximation for 
the genuine one; we can achieve an X rotation by Hill 
and Goan's scheme. The exact solution can be repre- 
sented by the eigenvalues and the eigenvectors of Eq. (jT3|l , 
which are given in AppendixlBl Adopting the typical 
value of the transverse magnetic field, B ac = 2.5 mT, 
and choosing Ox = 7r/4, we calculate F(\4>)i) for three 
values of the HY: A/A = 0.75, 0.5, 0.25, and differ- 
ent two initial states. The X rotation for Ox = tt/4 
plays an essential role in the Hadamard gate Ur- Uh = 

_ ie ™ I /4 e »a I /4 f ,™ z /4 > We show the regults m Xa _ 

ble [H] They imply that the values of F(\(f>)i) are al- 
most on the order of 10 1 , regardless of the initial states 
and the magnitude of the HY; the X rotation given 



TABLE II: Calculated values of the fidelity with 8 X = 
tt/4 and B ac = 2.5 mT. 







nm 


A/A 


\<P)i = |0}< 


\<P)i = (|0> ( + \l)i)/V2 


0.75 


0.72514 


0.70458 


0.5 


0.72458 


0.70454 


0.25 


0.73390 


0.70526 



by Eq. i|16|) involves a serious error, whose probability, 
1 — max|,M,. [F(|(/>)i)1, is estimated at about 10 _1 . Ac- 
cording to Ref. |15| , the error probability in terms of the 
fidelity (i.e., 1 — max^) . [-F(|</>)i)]) is estimated about at 
10~ 5 in the cnot gate. We will need to discuss the pos- 
sibility of quantum error correction [2j| . However, it is 
certain that the error in the X rotation in Ref. 15] is 
vastly larger than the one estimated there. 

The presence of the error means that F[ % mix in Eq. ijT3|) 
shouldn't be regarded as a perturbation; the bases where 
H d is a diagonal matrix are not suitable for the analysis 
of the Schrodinger equation (|12f> . The above considera- 
tion is self-evident, investigating the time evolution op- 
erator under the Lamour resonance condition i|14|) . The 
dynamics by H d , which is relevant to the X rotation, is 
characterized by the energy scale vgB ac . However, we 
usually find vg < fig. Consequently, the mixing effect by 
Hmix becomes dominant. We also recognize the domi- 
nant effect of H l mix from the following viewpoint of esti- 
mating errors. If we simply apply the ordinal perturba- 
tion method to the calculation of the eigenvectors of H l rot 
in Eq. I|17fl . we obtain the following result for \<fi)i = |0); 
under the condition l|14|) : 



mow = 




' R 2 

O I — 

1 B 2 



Here, we regard H d as the free part in (|13|) . We define 
the difference between the lowest eigenvalue of H d and 



the one of HI 



as Aq. 



We also define the difference 



between the second lowest eigenvalue of H d and the one 
of H l rot as Ai. Indeed, we find t x (A[ - A l )/fi « 10~ 6 
for Aj Aq = 0.5 and Ox — tt/4. Then, we approximately 
get F(\0)i) « 1 for A/A = 0.5 and X = tt/4; this result 
is clearly different from the rigorous one. 

In summary, we have shown that the X rotation in 
Ref. 0| includes a considerable error, which is generated 
by H % mix . The mixing Hamiltonian H l mix never appears 
under the assumption that the electron spin state is al- 
ways down, and this term should not be regarded as a 
perturbation for H l d . 



C. The controlled— Z gates 



In SecJTHBl we explained that there is considerable 
error in the X rotation carried out according to Hill and 



7 



Goan's scheme. Generally, spin-flip operations play an 
important role in controlled operations. In Ref. |15|, they 
are frequently used even in the controlled-Z gates. Here, 
we show that the error in the controlled-Z gate mainly 
occurs in the spin-flip operations. 

First, by means of the adiabatic theorem, we calcu- 
late the time evolution operator for two qubits during an 
interval tc, 



U c {tc) = T 



exp 



H c (t)dt 



(18) 



keeping Ai(t) = A 2 (t) = A n , turning off the transverse 
magnetic field, and adiabatically varying J(t). The total 
operation time of J(t) is defined as tc- In order to get 
an analytical expression of the operator, we take the fol- 
lowing profile of J(t), introducing four parameters: The 
maximum value of J(i), J c , the period th when J(t) keeps 
J c , the former and the latter periods when J(t) varies 
adiabatically, t a (i.e., tc = t a + U L ), and the smoothness 
of profile touching to the constant line, r'. The assumed 
time dependence of J(t) is as follows: 



J(t) = { Jc 

Jcfj(T C 



(19) 



(0<r<i) 

( T C - 2 - T - T Cj 

where r is a scaled time r = t/t ai tc = tc/t a , and 
Th = th/t a (see Fig.|3J). The function /j(r) in Eq. (JUJ is 
defined as, 



fj(r) = e(r' -t)-t 2 + e(r-r') 



(l"2r) 2 



1 - 2r' 



The initial state and final state for two qubits will be 
spanned by \v\), \v + ), and \vi), because J(0) = 

J(t c ) = and Ai(t) = A 2 {t) = (0 < t < t c ). Then, 
the time-evolution operator 1(18(1 is given by 



U C (U 



e- l ^ + ^\v k ){v k \+R c . (20) 



fc=l, +, -,4 



The phase a k and f3 k in Eq. ((20(1 are related to the pro- 
cesses in which J(t) varies (i.e., < r < 1/2 and 
tc — 1/2 < r < tc) and in which J(t) = J c (i.e., 
1/2 < r < tc — 1/2), respectively. The analytical ex- 
pressions for them are shown in Appendix E] The last 
term in Eq. 1(20|) . like the third term in Eq. |JS}, repre- 
sents the deviation from the adiabatic approximation and 
mixes the states relevant to quantum computation with 
the irrelevant states. This gives a dominant contribution 
near the level crossing region around J c = /is-B/2 (see, 
Fig.^l, where we have to apply carefully the adiabatic 
approximation. In this paper, we investigate Eq. (|18fl in 
the regime where the level crossing is not included; we 
can safely neglect Rc in Eq. (|20[1 . Actually, we confirm 
that the adiabatic approximation works very well, solving 
numerically the time-dependent Schrodinger equation. 

Secondly, we explain how to connect Uc{tc) with a 
controlled-Z gate. Hereafter, we assume that the spin- 
flip operations are realized by using Eq. 1(16(1. First of 



0.5 - 




FIG. 3: Profile of J(t) (r h = t h /t a = 2.0 and r' = 0.2). 



all, we rewrite Eq. as Uc(tc) — U a dU s t, where 
Uad = Efee- Mfc k)KI and U st = £ fe e^* k)<^|. 
Next, we investigate U s t in more detail. Introducing 
1, = |0)i,i(0| + |l)i,i(l|, a iz = |0)i,i(0|-|l)i,i(l|, a ix = 
|0)i,i{l| + |l)i,i(0|, ando-iv = -»(|0)£,i(l| -11)^(01), we 
find that U s t is given by 



where V = e~ lc e'^" 1 c = (/3i+/3+ + /3_+/3 4 )/4, 

S s = (/3i-/3 4 )/4, 8c = (/3+-/3_)/4, and S' c = {fit- 13+ - 
13- + AO/4. Note that this operator contains spin-flip 
operations. Here, we define the unitary operator W as 

W = (U H ®U H )V^U st (U H ®U H ) 

= e~ **C(01 z®Ol Z+&1 yt$&2 y) e -iS' c <Jl x<g>CT2 x 

where Uh is the Hadamard gate. Then, according to the 
idea shown in Ref. , we find that the controlled-Z gate 
UcziPcz) of the angle CZ is achieved as follows: 



= e 

= W(±i ® a 2z )W(t 1 ® a 2z ) 

x(l l( g)e 4 ^ CT2 ^. (21) 

Note that the several basic identities (Ji Z ai X ai Z = —<Ji X , 



'i z u i y t z 



-a ly , and [cti a <g> a 2 a , ai b ®a 2b \ = (a, b 



x, y, z) play an essential role in the derivation of Eq. 1)21(1. 
The angle 9 CZ is given by 8 CZ = 4S C = 0+ — from 
which th is given by, 



th = 



2 J c + y/(e - 2 J c ) 2 + AA 2 Q - + 4^ 



(22) 



Here, let us recall that the physical time-evolution is 
given by Uc{tc)- The unitary operator W in Eq. l|2*T|l 
is rewritten in the term of Uc(tc) as follows: 

W=(U H ® U H )V^ul d U c {tc){U H <g> U H ). (23) 

The unitary operator V is composed of several Z rota- 
tions. 

Now that we find the connection between Uc(tc) and 
the controlled-Z gate, we should discuss the effect of U a d- 
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Generally, the unitary operator U a d should include two 
qubit operations. We would like to represent U ad as sin- 
gle qubit operations, in particular, Z rotations, because 
they have been constructed in Sec. lIII Al and their errors 
are very small. We find that U a d is equal to the single 



qubit operation, (Ji z ®&2 z, if the following conditions are 
fulfilled: a + — a± = tt + 2to + 7t, a_ — (Xi = tt + 2to_7t, 
and «4 — ai — 2m^ix (m+, m_, TO4 S 1 and TO4 7^ 0). 
These three conditions are written by 



-v4 2 + 4Ag-e + 


2g nf i n B + 2 j} 


>)-^ +, (r) 


dr 


1 + 2m + 


-2e + Ag n p n B 


+ 2A + 2J? 




r)-4°' + \r) 


dr 


2rriA L 



+ 2^5 - 2 /* { V[ e -2J(r)]2+4A§ + J(r) + +) (r)} 



dr 



ta = 



-2e + 4 ffnM „S + 2A + 2 f Q 5 [j(r) - e{ ' +) (t) 

2TO47r/l 



1 + 2m_ 
2m4 



-2e + 4g n p n B + 2A + 2 f* J(r) - e[° ' +> (r) 



dT 



(24) 

(25) 
(26) 



The analytical expression for E{ ' ' (r) is shown in Ap- 
pendix ^ We find that Eqs. JMJl and J55j include two 
unknown parameters J c and r', besides the free param- 
eters m+, and 7724. We should choose suitable val- 
ues for the free parameters that bring about a satisfac- 
tory solution. First, we numerically find that both the 
right-hand sides of Eqs. Ij24(l and (|25J) are close to 1/2 
for < J c < e/2 and < t' < 1/2. Thus, we choose 
m + = rn_ = and 7714 = 1. After this consideration, we 
numerically solve the simultaneous systems of equation 
composed of Eqs. it2"4"|) and lj2"5|) by the Newton-Raphson 
method. Then, we obtain t a by substituting the value 
of ?TJ4, J c , and t' for Eq. 126f) . Finally, the remaining 
parameter th is determined by Eq. Q22J) . 

We show the typical solutions for 8 CZ — tt in Table ITTll 
The case of 9 CZ — n is very important, because a cnot 
gate is essentially composed of U cz (tt) and Hadamard 
gates. We find that if we choose a large value for J c , 
the corresponding values of t' and th decrease; the op- 
eration time of the controlled-Z gate becomes short. On 
the other hand, the value of t a for each J c is almost the 
same and greatly shorter than th ; it suggests that we can 
safely regard the adiabatic controlling part of the profile 
of J(t) (i.e., < r < 1/2 and r c - 1/2 < r < r c ) as an 
instantaneous process within the total controlling one. 

Let us summarize the above discussion. First, we ob- 
tained Eq. (|21(l , which represents the relation between 
the logical controlled-Z gate and the physical time evo- 
lution, assuming the validity of the adiabatic approxi- 
mation and the realization of the spin-flip operations by 
Eq. I|21ll . The result is essentially equivalent to that of 
Ref. 1 1 5| . except for the presence of U a d- Here, we have 
shown that U a d is equivalent to the single qubit opera- 
tion, <7i z ® (T2 z, if we choose a suitable set of the pa- 
rameters (J c , t', t a , th). We find that one of the two 
assumptions, the adiabatic approximation, is valid, be- 



TABLE III: Calculated values of J c , r', t a , and th with 9 CZ — 
tt, m+ = 0, m_ = 0, and 1714 — 1. 



J c /e 


/ 

T 


t a (ns) 


th (ps) 


0.1003 


0.1085 


5.391 


33.1 


0.1988 


0.0916 


5.391 


12.59 


0.01 


0.2203 


5.392 


40.44 



cause we choose a value of J c which is far from the level 
crossing point. Unfortunately, the spin-flip operations 
by Eq. I|21l) don't quite work. The main origin of the er- 
ror in the controlled-Z gates lies in the part relevant to 
the spin-flip operations (e.g., the Hadamard gate). 



IV. SUMMARY 

We have re-investigated the operating schemes of 
quantum gates associated with universal quantum gates 
(i.e., Z rotations, X rotations, and controlled-Z gates), 
for Kane's model in detail. We chose the suitable 
computational bases out of the eigenstates of H l (t) 
with fixed time. Our choice is different from those in 
Refs. 0, 0, 0, in which the electron spin state is as- 
sumed to be always down. The bases in the latter case 
are fragile for the time-evolution due to H l (t). The 
phase-shift operations (i.e., Z rotations) for the i th qubit 
are shown to be constructed with an extremely low er- 
ror probability, adiabatically varying the value of Ai(t) 
with the fixed other parameters: J(t) = 0, Aj(t) = Aq, 
and B ac = 0. The most important result is that a 
considerable error for the X rotation in Ref. exists. 
The physical origin of such an error is the mixing ef- 
fect between the computational bases and the irrelevant 
states by B % mix , which is neglected in Ref. . This issue 
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emerged from discussion of Eq. i|12|) , not by perturbation 
but by a rigorous analysis. Finally, the controlled opera- 
tions in Ref. will not quite work, because the spin-flip 
operations are used in them. However, it is meaningful 
to check whether they contain other errors. In particu- 
lar, we have investigated the controllcd-Z gate, as it is 
an essential part of the cnot gate. We have shown that 
the main origin of the error in the controlled- Z gate lies 
in the parts relevant to the spin-flip operations. Accord- 
ingly, we will have to discuss recovery of the errors caused 
by the spin-flip operations. 

It is necessary for quantum computation to control 
quantum systems. In doing this, it is inevitable that 
some errors occur in the controlling processes, even if we 
do them carefully. Therefore, it is crucial to make clear 
the physical origin of the errors and to investigate meth- 
ods of correcting them. Kane's model is an attractive 
proposal for quantum computers. In this paper, we have 
presented the robust quantum computational bases for 
H % (t) time-evolution, but they do not quite work for the 
spin-flip operations. We will need to find more suitable 
computational bases for the operations of quantum gates. 
In addition, we will have to take account of the modified 
versions of Kane's model 0i E3- 



Acknowledgments 

The authors acknowledge H. Nakazato, K. Yuasa, T. 
Watanabe, and T. Shinada for valuable discussions. This 
work is supported partially by a Grant-in-Aid for the 
COE Research at Waseda University and that for Pri- 
ority Area B (No. 763), MEXT, by a Grant for the 21st 
Century COE Program at Waseda University, and by a 
Waseda University Grant for Special Research Projects 
(No.2004B-872). The authors thank the Yukawa In- 
stitute for Theoretical Physics at Kyoto University. 
Discussions during the YITP workshop YITP-W-04- 
14 on "Chaos and Nonlinear Dynamics in Quantum- 
Mechanical and Macroscopic Systems" were useful to the 
completion of this work. 



APPENDIX A: THE INSTANTANEOUS 
EIGENVALUES OF THE TWO QUBIT 
HAMILTONIAN 

We show the explicit expressions for the instantaneous 
eigenvalues of Hc(t) at the time t if A\(t) = A 2 (t)(= 
A(t)). The Hamiltonian Hc(t) is a 16 x 16 matrix and 
is decomposed into several submatrices by the system's 
symmetry (see Eq. (|2J)). Hereafter, we abbreviate the ar- 
gument t. We only write the eigenvalues relevant to the 



computation bases: 

-\h~ ei + ^ <A1) 
4°'" = -A-\-Jt^ 

< A2 » 

E[-^- 1] = -J-e + 2g ni i n B 

- v / (e-2J) 2 + 4A 2 , (A3) 

= -J-e + 2g n fi n B 

+ V(e-2J) 2 + 4A 2 , (A4) 

Ej~ M) - J-e + 2g nf i n B- V / e 2 +4 J 4 2 , (A5) 

E { ~ 2 ' 1] = J - 2e + 4g n (j, n + 2A. (A6) 

We define £ in Eqs. (|A1|) and (|A2|) as a real root of the 
following cubic equation: 

z 3 - az 2 - 4c 3 z + 4c 3 ci - c\ = 0. (A7) 

Taking account of Cardano's formula 26], we easily ob- 
tain a real root of Eq. I|A7J| . The coefficients in Eqs. I|A1(I . 
(|A2|) . and (|A7I) . Ci, c 2 , and c 3 , are given by, 

Cl = -6J 2 + 4JA-22A 2 -4e 2 , 

c 2 = 8J 3 -8AJ 2 + 8(A 2 - e 2 )J + 8e 2 A + 2AA 3 , 

c 3 = -3J 4 + 4A/ 3 + {UA 2 + 12e 2 ) J 2 

+ (-6(L4 3 + 8e 2 A)J + l2e 2 A 2 + 45A 4 . 

Using Eqs. fTi^l and l|Al|l - ljA6|l . the phases ak and (3k in 
Ea. (l20|) are given as follows: 

«l = ^^ 1/2 ^°' 1) (J(r))dr, 

o± = ^^ 1/2 ^- 1 ' ±1) (J(r))dr, 

a 4 = ^J^ 2 E[- 2 ^(J( T ))dr, 

01 = f^'Vc), 

A = \E^\J C ). 

APPENDIX B: THE SOLUTION OF THE 
EIGENVALUE PROBLEM FOR THE SINGLE 
QUBIT HAMILTONIAN IN THE ROTATING 
FRAME 

Let us consider the solution for the eigenvalue problem, 
Hl ot \w k )i = hfll\w k )i (k = 0, 1, 2, 3). The eigenvalues 



10 



are given by 

mi 



m\ 



tei\ 



:Vv - di 



--^Vv-di 



1 , 2d 2 



-s/v-di 



1 , 2d 2 



m\ = -Vv-di 



1 , 2d 2 



(Bl) 



(B2) 



(B3) 



(B4) 



We define 77 in Eqs. (|B1|) - (|B4(1 as a real root of the fol- 
lowing cubic equation: 



diz 2 - Ad 3 z + Adxd 3 



d\ 



0. 



(B5) 



We easily find a real root of Eq. (|B5J| through Cardano's 
formula |2q . as in Appendix 1X1 The coefficients, d\ % d 2 , 
and d 3 in Eqs. (|Bl(l - i|B5(l . are given by 



di 
d 2 
d 3 



A) ^2 



(E l + E l 2 ) 2 + (El - hio ac ){El + hw ac ) 
-2(fi B B ac ) 2 -2{g n ^ n B ac ) 2 , 
(E l + E\) [E l Ei - (E{ - hw ac )(El + hw ac )} 

EiEi{E\ - hw ac )(E* 3 + hw ac ) + ( M | - g 2 y n ) 2 B 2 ac 
-{E\-hw ac ){Ey_ e +E^ 2 _ e )B 2 lc 
-(El + hj ac )(E^ 2 +E^ 2 )B 2 ac . 



The eigenvector \wk)i associated with the eigenvalue Ml\ 
is given by 



\Wk)i 



v_ 6 El - mi E\ - m\ 



lis vgB a 



HeBac 



E\ + hw ac - tel l k E o 



u (A))i 

mi 



E^ tel ac — ^^l E\ 



mi 



veB a 



H-eB a 



HeB a 



v-e 



E l 3 + hw ac 



m% e\ - m\ El 



tel\ 



\MA))i 



El - tell 



H-e E\ - tell 



^eB a 



H—eB ac 



\MA)) l 



where Af is a normalization factor. 
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